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I.  INTRODUCTION 


The  passive  Rb  gas  cell  atomic  clock  is  presently  the  most  often  used 
atomic  frequency  standard  in  high-performance  systems  because  of  its  par¬ 
ticular  advantages  of  small  size,  low  weight,  low  power  consumption,  and 
excellent  short-term  performance.  In  this  standard  the  Rb  atomic  resonance 
frequency  associated  with  the  hyperfine  ground  state  splitting  is  used  as  a 
reference  to  control  the  frequency  of  a  quartz  crystal  oscillator  (see  Fig.  1 
and  Ref.  1).  Servo  control  circuitry  is  used  to  transfer  the  stability  of  the 
microwave  atomic  transition  to  the  crystal  oscillator,  and  for  times  greater 
than  the  time  constant  of  the  feedback  loop,  the  crystal  oscillator  frequency 
stability  reflects  that  of  the  atomic  system.  The  Rb  hyperfine  resonance  is 
optically  detected  using  the  standard  technique  of  optical  pumping  double 
resonance,  where  the  source  of  optical  pumping  radiation  is  typically  a  Rb®^ 
discharge  lamp.  For  effective  optical  pumping,  the  output  of  the  discharge 
lamp  is  filtered  by  transmission  through  a  vapor  of  Rb^  atoms.  The  filtering 
may  be  carried  out  using  a  discrete  separated  filter  cell  (SFC),  as  is  shown 
in  Fig.  1,  or  by  combining  the  filter  and  absorption  cells  into  an  Integrated 
filter  cell  (IFC).1'2  In  either  case  the  atomic  microwave  signal  appears  at  a 
photodiode  as  a  small  light  intensity  modulation  component  resting  on  a  large 
DC  background.  In  this  report,  a  nonempirical  model  of  the  gas  cell  frequency 
standard  is  developed,  with  the  aim  of  putting  the  operation  and  performance 
of  the  gas  cell  clock  on  a  sound  theoretical  footing. 

In  past  years,  a  number  of  other  authors  have  investigated  Rb  clock  per- 

1  3-7 

formance  by  developing  clock  models  of  varying  complexities.  *  We  would  be 
remiss  if  we  did  not  acknowledge  the  foundation  these  studies  have  provided 
for  our  own  work.  The  present  model,  though,  has  a  number  of  significant 
differences  from  the  previous  investigations.  All  previous  models  were 
limited  to  short-term  performance  analysis,  using  Allan  variance  averaging 
times  of  less  than  lO^sec.  During  this  time  regime,  shot  noise  at  the 
photodetector  limits  clock  performance.  The  current  analysis  is  the  first 
that  is  valid  for  times  greater  than  lO^sec.  This  was  accomplished  by 
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Fig.  1.  (a)  Schematic  Diagram  of  a  Passive  Rubidium  Gas  Cell  Frequency 

Standard;  (b)  Energy  Level  Diagram  of  Rb®  Exhibiting  the  Radiative 
Transitions  of  Interest 


explicitly  including  the  effects  of  discharge  lamp  intensity  fluctuations, 
which  are  transferred  into  output  frequency  variations  via  the  light-shift 
effect.  Additionally,  the  model  is  the  first  which  incorporates  the  influence 
of  an  optically  thick  Rb  vapor  along  with  the  diffusion  of  optically  pumped 
atoms  to  the  walls  of  the  absorption  cell  into  the  calculation  of  clock 
performance.  Inclusion  of  these  processes  very  dramatically  affects  the  axial 
Rb  hyperfine  polarization  within  the  cell  and  how  each  segment  of  the  cell 
contributes  to  the  observed  hyperfine  resonance  signal.  With  the  advent  of 

O 

diode  laser  sources  for  optical  pumping,  practical  gas  cell  clocks  based  on 
alkalies  other  than  Rb,  e.g.,  Cs,  may  be  envisioned.  Consequently,  the 
current  analysis  has  been  carried  out  for  arbitrary  nuclear  spin,  permitting 
analysis  of  alternative  gas  cell  clocks. 

This  report  has  been  divided  into  several  sections.  In  the  analysis 
section,  the  feedback  analysis  of  the  servo-control  loop,  along  with  the  phy¬ 
sical  analysis  and  computational  procedures  yielding  the  Rb  hyperfine 
lineshape,  are  developed.  Integration  of  these  components  yields  the 
functioning  model.  In  the  model  verification  section,  the  validity  of  the 
model  is  established  by  modeling  the  performance  of  a  prototype  "state-of-the- 
art"  Rb  standard  »  u  manufactured  by  EG&G,  Electronic  Components  Division,  for 
use  on  NAVSTAR/ global  positioning  system  satellites.  In  the  final  section, 
implications  of  this  clock  model  are  discussed  along  with  potential  avenues 
for  Rb  clock  performance  Improvement. 


II.  ANALYSIS 


A.  FEEDBACK  ANALYSIS  OF  SERVO  CONTROL  CIRCUITRY 


The  function  of  the  servo  control  circuitry  is  to  transfer  the  stability 

of  the  microwave  atomic  transition  to  the  crystal  oscillator.  The  feedback 

analysis  of  this  system  provides  a  means  of  understanding  how  the  electronics 

and  the  physics  package  affects  clock  performance  at  all  analysis  times.  In 

Fig.  2  the  clock  is  displayed  in  a  functional  block  diagram  to  which  standard 

linear  feedback  analysis  may  be  applied.  The  analysis  is  carried  out  in  the 

frequency  domain,  and  the  inputs  are  the  power  spectral  densities  of  the 

various  system  noise  sources.  Scr(f)  and  Sph(f)  are  the  spectral  densities  of 

y  y 

the  fractional  frequency  fluctuations  of  the  crystal  oscillator  and  the 
physics  package,  respectively.^  Other  authors  have  carried  out  most  of  the 
required  feedback  analysis  previously . ^ ^ ^  The  spectral  density  of  the 
fractional  frequency  fluctuations  of  the  entire  clock,  voltage  controlled 
quartz  crystal  oscillator  (VCXO)  locked  to  the  atomic  system,  SRb(f),  is  given 
by 


S*b(f) 


1  +  (f /fn): 


Syr(f)  ♦ 


Sph(f) 


i  +  Cf/fny 
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where  R  is  the  effective  impedance  of  the  current  to  the  voltage  converter/ 
amplifier,  M  is  the  frequency  multiplier's  multiplication  factor,  and  m  is  the 
atomic  system  response  relating  a  microwave  frequency  detuning  Af  to  a 
detected  error  signal  current  i, 


Effectively,  o  is  found  to  be  the  slope  of  the  discriminator  pattern,  propor¬ 
tional  to  the  second  derivative  of  the  Rb  hyperfine  llneshape  evaluated  at  the 
Rb  resonance's  central  frequency.  This  parameter  is  obtained  through  modeling 
of  the  physics  package  (see  subsections  II. B  and  II. C).  The  VCXO  response 
li  relates  a  change  in  control  voltage  AV  to  a  change  in  the  crystal's  output 
frequency  Af^  , 


Af 


VCXO 


WAV 


in 


(3) 


In  Eq.  (lb),  tj  is  the  time  constant  of  the  integrating  filter  which  has  a 
transfer  function  of  the  form  l/^wiftj).  Equation  (la)  indicates  that  the 
spectral  density  of  the  clock's  fractional  frequency  fluctuations  for 
frequencies  above  fn  is  determined  by  the  crystal  oscillator  and  for 
frequencies  below  fQ  by  the  physics  package.  In  the  time  domain,  this  implies 
that  for  times  less  than  the  loop  attack  time  (1/2*^),  the  clock  frequency 
stability  is  specified  by  the  crystal  oscillator;  for  times  greater  than  the 
attack  time,  frequency  stability  is  defined  by  the  atomic  system  (physics 
package ) . 

The  principal  difference  between  previous  analyses  and  the  present  one  is 
the  explicit  inclusion  of  the  effects  of  the  lamp-induced  light  shift1^’1^  of 
the  Rb  hyperfine  transition  frequency.  It  is  clear  that  lamp  intensity 
fluctuations  at  frequencies  less  than  fn  are  mapped  directly  into  variations 
of  the  clock  output  frequency  via  the  light-shift  effect.  Recent  experimental 
studies  have  shown  that  for  the  EG&G  Rb  frequency  standard,  these  light- 
induced  frequency  fluctuations  specify  the  long-term  frequency  stability  of 
the  standard.^  Further,  these  lamp-induced  frequency  fluctuations  should  set 
a  fundamental  limitation  to  the  long-term  frequency  stability  of  Rb  clocks  in 
general.  Including  the  light-shift  effect,  the  physics  package's  spectral 
density  of  fractional  frequency  fluctuations  takes  the  form 


sjh(f)  -  [kJs  sj(f)  + 


COr 


where  Che  light-shift  coefficient  KLS  relates  a  change  in  clock  fre¬ 
quency  Af^,,  and  the  nominal  clock  frequency  f°^  to  a  change  in  light 
Intensity  AI,  and  the  nominal  light  intensity  I 


4£LS  /fRb  ■  <4I'V 


In  Eq.  (4),  S*  (f)  is  the  power  spectral  density  of  the  fractional  light 
intensity  (AI/I^)  fluctuations.  The  second  term  on  the  right-hand  side  of  the 
equation  is  the  result  of  shot  noise  at  the  photodetector,  where  e  is  the 
electronic  charge,  iQ  is  the  total  lamp-induced  photocurrent,  and  f®  is  the 
Rb  hyperfine  resonance  frequency. 

Modeling  of  the  physics  package  supplies  iQ,  while  S*(f)  can  be 
experimentally  determined  from  lamp  intensity  fluctuations  recorded  as  a 
function  of  time.  Whereas  In  principle  the  light-shift  coefficient  can  be 
calculated  from  fundamental  parameters,  in  the  present  analysis  an  experi¬ 
mentally  measured  value  is  used.  Equation  (4)  explicitly  Indicates  how  lamp 
Intensity  fluctuations  and  the  light  shift  can  impact  clock  performance.  The 
Importance  of  minimizing  the  light-shift  coefficient  as  well  as  the  design  of 
lamps  yielding  extremely  stable  Intensities  are  apparent. 

The  spectral  density  of  fractional  frequency  fluctuations  is  related  to 
the  Allan  variance  o^(t)  via  the  following  integral  transformation:11 


(?) 


2  /  S*\f) 


sin  (irff) 
(*rf )2 


where  t  is  the  Allan  variance  frequency  averaging  time.  For  the  EG&G  Rb  fre¬ 
quency  standard,  S*(f)  was  of  the  form16 


(6) 


S*(f)  -  A  +  K 
y  f 


which  Is  the  sum  of  white  and  random  walk  contributions.  Crystal  oscillators 
display  fractional  frequency  spectral  densities  of  the  form  ^ 


s“(f)  -  j  +  Df  +  Ef2 


(7) 


Cbrabining  this  information  with  Eq.  (2),  we  observe  that  the  following 
integrals  must  be  evaluated  when  obtaining  the  Allan  variance  over  all 
averaging  times: 


—  f  ®  A 

A(t)  -  /  - =■  sin  (rtf)  df  (8) 

0  1  +  (f/fnr 

a  -  -  4,  -2,  -1,  1,  2 


To  permit  rapid  evaluation  of  Allan  variances,  particularly  for  averaging 
times  near  the  loop  attack  time,  these  Integrals  were  evaluated  analytically 
using  complex  contours.  The  results  of  these  evaluations  are  tabulated  in  the 
appendix. 

With  the  evaluation  of  the  integrals,  the  Allan  variance  may  be  readily 
extracted  from  the  modeled  clock's  spectral  density  of  fractional  frequency 
fluctuations.  This  spectral  density  is  a  function  of  crystal  and  loop 
parameters  as  well  as  inputs  describing  the  performance  of  the  physics 
package.  In  the  following  paragraphs,  the  analysis  applied  to  the  physics 
package,  yielding  the  total  detected  photocurrent  and  the  atomic  system 
response,  m,  is  described. 
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B 


ATOMIC  SYSTEM  ANALYSIS 


Figure  la  exhibits  the  basic  physics  package  elements  of  an  optically 

pumped  gas  cell  atomic  frequency  standard:  the  lamp,  the  filter  cell,  the 

resonance  cell,  and  the  photodetector*  In  combination,  the  lamp  and  filter 

cell  produce  a  spectral  emission  that  preferentially  excites  atoms  In  the 

2 

resonance  cell  out  of  one  of  the  5  ground  state  hyperflne  multlplets. 

For  example,  the  optical  emission  lines  corresponding  to  transitions  from 

F  »  2  (i.e.,  transition  numbers  3  and  4  in  Fig.  lb)  may  be  effectively  blocked 

by  the  filter  cell,  so  that  only  those  atoms  in  the  F  ■  1  hyperflne  multlplet 

can  absorb  light.  Because  an  excited  atom  can  decay  from  the  5 ^3/2  or  5  Pj^ 

state  Into  either  ground  state  hyperflne  multlplet,  the  successive  action  of 

several  optical  absorptions  and  reemlsslons  will  lead  to  a  transfer  of 

population  from  one  hyperflne  multlplet  into  the  other.  Furthermore,  because 

this  process  of  optical  pumping  depopulates  the  optically  absorbing  hyperflne 

multlplet,  the  equlllbriun  light  Intensity  transmitted  by  the  resonance  cell 

and  detected  by  the  photodiode  will  be  at  a  maximum.  If  microwaves  of  the 

appropriate  frequency,  6835  MHz,  are  now  applied  to  the  cavity,  atoms  in  the 

overpopulated  hyperflne  multlplet  will  be  Induced  into  returning  to  the 

optically  absorbing  hyperflne  multlplet.  This  appropriate  frequency  Is  then 

detected  as  a  decrease  in  the  light  intensity  transmitted  by  the  resonance 

7  8 

cell.  Because  the  Q  of  the  atomic  Interaction  with  the  microwaves  is  10  -10  , 
the  decrease  in  transmitted  light  Intensity  can  be  used  to  discriminate 
against  very  small  changes  In  the  6835  MHz  microwave  frequency.* 

From  the  outline  of  gas  cell  frequency  standard  operation  presented 
above.  It  Is  clear  that  an  accurate  calculation  of  the  response  of  the  trans¬ 
mitted  light  intensity  to  variations  in  the  microwave  frequency  is  a  primary 
concern  for  a  model  of  clock  performance.  Typically,  this  problem  Is 
approached  in  one  of  three  directions.  Quite  often  one  attempts  to  describe 
the  atomic  system  as  a  two-level  atom  interacting  with  the  microwaves,  an 
approach  that  has  the  advantage  of  being  easy  to  work  with  and  Intuitive.  In 
particular,  when  one  is  concerned  with  obtaining  insight  into  difficult 

problems,  for  example,  the  Influence  of  microwave  modulation  on  the  atomic 
19 

llneshape,  a  two-level  atom  description  is  the  most  reasonable  first  step. 


However,  this  approach  is  quantitatively  incorrect  and  can  even  be  in  quali¬ 
tative  disagreement  with  reality  under  normal  frequency  standard  operating 
conditions.  An  alternative  approach  is  to  start  from  first  principles  and, 
as  much  as  possible,  to  rigorously  calculate  collisional  interactions,  diffusion 
effects,  and  optical  transmission  profiles  through  the  resonance  cell.^-^ 

This  program  will  yield  quantitatively  correct  results;  however,  in  light  of 
the  number  of  interrelated  physical  processes  that  occur  in  the  gas  cell 
frequency  standard,  it  appears  to  be  an  unreasonably  complex  procedure  for 
predicting  clock  performance.  Furthermore,  as  more  physical  processes  are 
incorporated  into  the  model,  this  rigorous  approach  will  soon  lose  the 
intuitive  nature  that  makes  the  two-level  atom  so  attractive.  Obviously,  a 
middle  approach  is  desirable,  where  important  physical  processes  are 
recognized  and  treated  fairly  rigorously;  whereas  less  important  processes  are 
treated  by  suitable  simplifying  approximations.  In  this  regard,  the  starting 
point  of  the  model  to  be  presented  subsequently  is  the  generalized  Vanler 
theory  of  the  0-0  hyperfine  lineshape  in  optically  pumped  alkali-metal 
vapors. 6*20,25 

In  the  generalized  Vanier  theory,  all  the  Zeeman  sublevels  of  an  alkali's 
ground  state  are  considered.  However,  assume  that:  (1)  all  Zeeman  sublevels 
of  a  given  hyperfine  multlplet  interact  with  the  light  equivalently;  (2)  all 
Zeeman  sublevels  of  the  ground  state  have  an  equal  probability  of  being 
populated  as  a  result  of  decay  from  either  the  or  excited 

states;  and  (3)  the  microwaves  only  induce  transitions  between  the  (F  •  2, 

*  0)  and  the  (1,0)  ground  state  Zeeman  sublevels.  Using  these  assumptions, 
the  density  matrix  equations  describing  the  evolution  of  the  ground  state 
Zeeman  sublevels  are  considerably  simplified.  Additionally,  relaxation  of  the 
ground  state  Zeeman  sublevels  is  included  phenomenologically  into  the  model  by 
considering  "dark"  longitudinal  and  transverse  relaxation  rates  y^  and  y^, 
respectively.  Thus,  one  does  not  complicate  the  physical  picture  by  differen¬ 
tiating  among  uniform,  electron  randomization  and  Zeeman  relaxation 
mechanisms.  [At  the  present  time.,  there  appears  to  be  no  compelling  reason 
for  considering  relaxation  more  rigorously  (see  Ref.  20).]  In  the  remainder 
of  this  section,  we  will  use  the  generalized  Vanier  theory  to  construct  a 


model  for  Che  transmitted  light  intensity  as  a  function  of  microwave  frequency 
in  a  gas  cell  frequency  standard  that  is,  possibly,  optically  thick.  In  this 
regard,  our  procedure  for  constructing  the  model  will  be  to  consider  each 
Important  clock  element  and  physical  process  individually  before  incorporating 
it  into  the  evolving  clock  model. 

1.  LAMP  AND  FILTER  CELL 

We  assume  that  the  lamp  light  entering  the  resonance  cell  can  be 
decomposed  into  five  parts:  four  Gaussian  lineshapes  of  full  width  Av_ 
corresponding  to  transitions  1-4  in  Fig.  lb,  and  a  component  corresponding  to 
the  lamp's  buffer  gas  light,  which  contributes  to  the  shot  noise  at  the  photo¬ 
detector.  Defining  e  as  the  ratio  of  total  buffer  gas  light “to  total  rubidium 
light,  we  have  for  the  optical  spectrum  entering  the  resonance  cell 

4  ...  41n2( v  -  v. )2 

I(v)  "  J  ‘  VBg)  +  - 2  exp  ( - 2 - ^  (9) 

i-1  irAvG  Avc 


where  Ilo  is  the  total  entrance  intensity  of  line  i,  and  is  the  center 


*.U 

frequency  of  the  i  spectral  component.  The  presence  of  the  Dirac  delta 


function  A(v  -  v__)  is  a  simple  means  of  including  lamp  buffer  gas  light  into 

OVI 


the  analysis.  The  I^Q  are  determined  by  the  properties  of  the  lamp  and  the 
action  of  the  filter  cell,  both  of  which  have  recently  received  considerable 
study. ^“31  Our  model  should  be  applicable  to  both  the  SFC  and  IFC  Rb  clock 
designs.  This  flexibility  was  Included  by  allowing  for  an  isotopic  mixture  in 
the  resonance  cell  and  by  maintaining  the  I^Q  as  input  parameters  to  the  model 
(so  that  the  I1(J  can  represent  the  direct  spectral  emission  from  either  the 
lamp  or  the  lamp/filter  cell  combination). 


2.  TRANSMITTED  LIGHT 


To  take  into  account  the  effects  of  optical  pumping  on  the  light 
transmitted  by  the  resonance  cell,  we  write  the  Bouger-Lambert  law  for  the 
spectral  components  as 
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3Ii(A,z) 

3z 


(10a) 


-  -  nA(z)[Rb87]aiIi(z)  (1  -  1,2) 

{H  -  nA(z)]  iRb87 1  +  (yy)  [Rb85]}  <JjIj(z)  (j  -  3,4)  (10b) 

Equation  (10b)  reflects  the  fact  that  spectral  components  3  and  4  are  absorbed 

87  oc 

by  not  only  Rb  atoms  but  also  Rb°  atoms,  which  would  be  present  if  an  IFC 
design  were  employed.  In  these  expressions,  n,(z)  is  the  fractional 
population  in  the  F  *  1  hyperfine  multiple!  at  some  axial  position  z  in  the 
resonance  cell,  and  for  some  detuning  A  between  the  hyperfine  transition 
frequency  and  the  applied  microwave  field;  [Rb87]  and  [Rb85]  are  number 
densities;  and  the  are  optical  absorption  cross  sections: 

°u  “  2roV  /**n2/(Av2  +  (ID 

tH 

where  f  is  the  u  transition  oscillator  strength,  rQ  is  the  classical  elec¬ 
tron  radius,  and  Av^  is  the  Doppler  broadened  absorption  linewidth.  (Note 
that  for  an  isotoplcally  mixed  resonance  cell,  we  assume  negligible  optical 
pumping  of  the  Rb  Isotope.  Considering  the  low  light  levels  involved,  this 
is  probably  a  reasonable  first  order  approximation.) 

We  note  that  the  present  model  of  gas  cell  clock  absorption  is  basically 
one  dimensional;  that  is,  in  cylindrical  coordinates  we  consider  the  trans¬ 
mitted  intensity  at  a  specific  axial  position  z  to  be  equivalent  for  all  r 
and  0.  In  actuality,  however,  we  know  that  this  is  not  the  case  in  presently 
manufactured  designs,  because  the  microwave  field  strength  in  the  clock  cavity 
can  have  both  a  radial  and  angular  dependency.5^'55  Thus,  as  applied  to 
presently  available  gas  cell  frequency  standards,  the  present  model  must  be 
considered  as  a  first-order  approximation.  It  is  our  belief,  though,  that 
corrections  resulting  from  the  radial  and  angular  distribution  of  microwave 
field  strength  within  the  clock  cavity  are  within  the  uncertainties  generated 
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by  the  imperfect  knowledge  of  the  input  parameters  required  by  the  model. 
Furthermore,  it  is  likely  that  advanced  gas  cell  frequency  standards  will 
employ  dielectrically  loaded  cavities  whose  microwave  fields  have  no  radial  or 
angular  dependencies,8^  or  photodetectors  that  only  sample  transmitted  light 
intensities  at  specific  values  of  r  and  e.8^  For  these  advanced  frequency 
standards,  we  expect  our  model  to  be  somewhat  more  rigorous. 

Integrating  Eqs.  (10a)  and  (10b),  the  intensity  of  the  rubidium  light  at 
any  axial  position  z  is  given  by 


I1(A,z)  -  Ilo  exp  {-  nA(z)  [Rb87]  a±z)  (i  -  1,2) 


(12a) 


Ij(A,z)  -  exp  (-{[1  -  hA(z)l  [Rb87]  +  (yj) 

where  n.(z)  is  the  fractional  population  in  the 
A 

averaged  over  the  vapor  length  traversed  by  the 


[Rb85]}  at)  (j  -  3,4)  (12b) 

F  *  l  hyperfine  raultlplet 
light : 


Vz) 


z 

j  nfi(z)  dz 
o 


(13) 


[For  convenience,  at  z  -  l,  n.(i)  will  simply  be  denoted  as  n..]  Thus,  the 

A  A 

quantity  of  Interest  for  determining  clock  performance  becomes 


I(A,t) 


Ilo  (e  +  exp  {— [Rb87 ] 
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♦  I 

j-3 


e  +  exp(-  {(1 


-  n&)  [R&7  1  +  (yj)  [Rb85]}ojl  )  j. 


(14) 


This  equation  is  deceptively  simple  in  appearance.  The  transmitted  light 


intensity  is  just  the  sura  of  several  exponential  functions.  In  actuality, 
the  terras  appearing  in  the  exponents  greatly  complicate  the  situation, 
because  they  depend  on  the  lt(z).  Thus,  the  equation  for  I(A,i)  is  trans¬ 
cendental  and  requires  special  consideration.  Note  that  the  transcendental 
nature  of  Eq.  (14)  results  solely  from  the  fact  that  under  normal  clock 
operating  conditions,  the  resonance  cell’s  alkali  vapor  is  optically  thick  to 

the  rubidium  lamp  light.  For  an  optically  thin  vapor,  rT"  would  be  a  function 

& 

of  the  given  I  . 

MO 

3.  THE  FRACTIONAL  POPULATION  n 

A 

In  the  generalized  Vanier  theory,  the  fractional  population  in  the  lower 
F  *  b  •  I  -  1/2  hyperflne  multiples,  where  I  is  the  alkali's  nuclear  spin  and 
diffusion  is  considered  phenomenologically,  is  given  by  the  expression20 


g.U  +  y.) 
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_ •  .  f  b  1  f  2  "1  , 
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where 


and 


IB 
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ri.  -  « *  v  ii  *  <if->  <Hr,,rl 
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(A  ♦  Y  j )  (B  +  Yj)  (g^B  ♦  gfeA  ♦  gYj) 
g((A  ♦  B)/2 1 2  ♦  A(B  -  A)  ♦  Yj{g((A  +  B)/2  ♦  Yj)  +  gfiB  +  g^} 


(16) 


(17) 


(18) 


In  these  expressions,  the  subscripts  a  and  b  refer  to  the  F  quantum  numbers  of 
the  hyperflne  multlplets  (e.g.,  a  -  I  +  1/2),  and  uij  is  the  microwave  Rabi 
frequency.  The  g's  refer  to  the  various  degeneracies  in  the  alkali  ground 
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state:  g  is  the  total  ground  state  degeneracy  equal  to  2(21  +  1),  gg  is  the 
degeneracy  of  the  upper  F  -  a  hyperfine  multiplet  equal  to  2(1  +  l),  and  gb  is 
the  degeneracy  of  the  lower  hyperfine  multiplet  equal  to  21  (for  Rb®7  g  *  8, 
ga  -  5,  and  gb  -  3).  The  terms  A  and  B  are  defined  as  optical  photon 
absorption  rates  for  the  F  ■  2  and  F  =  l  hyperfine  raultiplets,  respectively. 
Because  of  the  variation  of  optical  pumping  light  intensity  with  axial 
position,  the  magnitudes  of  these  quantities  depend  on  z.  The  calculation  of 
A  and  B  is  carried  out  subsequently.  The  quantities  which  depend  on  axial 
position  A,  B,  and  result  in  also  displaying  an  axial  dependence. 
However,  because  of  the  inclusion  of  the  effects  of  atomic  axial  diffusion, 
the  quantity  of  interest  n  has  a  different  axial  dependence  than  n ' ,  though, 
as  shown  below,  there  is  a  relatively  simple  expression  relating  them. 

Probably  the  stickiest  issue  for  the  generalized  Vanier  theory  to  handle 

is  that  of  atomic  diffusion.  In  uncoated  resonance  cells,  after  an  alkali 

atom  hits  the  wall,  the  population  density  distribution  among  all  of  the 

atom's  ground  state  Zeeman  sublevels  is  randomized.35  Thus,  the  measure  of 

hyperfine  polarization  <t*£>  is  zero  at  the  resonance  cell  walls.  Rigorously, 
2 

a  term  DV  p  must  be  included  into  the  density  matrix  rate  equations  in  order 

to  describe  diffusion,  where  D  is  the  diffusion  coefficient  of  Rb  atoms  in  the 

36 

resonance  cell's  vapor.  However,  in  the  generalized  Vanier  theory,  diffu¬ 
sion  is  considered  as  a  phenomenological  relaxation  mechanism,  contributing  a 

term  y...  to  both  y  and  y  .  If  R  is  taken  as  a  measure  of  the  resonance 
dif  12 

cell's  dimensions,  typically,  the  resonance  cell's  radius,  then  y...  *  D/R  . 

air 

This  procedure  for  including  diffusion  into  the  determination  of  n,  is 

A 

actually  quite  reasonable,  because  if  no  other  relaxation  mechanisms  were 
present,  n^  would  relax  in  the  dark  (i.e.,  in  the  absence  of  optical  pumping 
light  and  microwaves)  with  a  rate  of  this  order  of  magnitude.37  However, 
diffusion  is  fundamentally  different  from  other  "bulk”  relaxation  mechanisms 
like  Rb-buffer  gas  and  Rb-Rb  collisions,  because,  in  steady  state,  diffusion 
leads  to  a  spatial  distribution  of  <!»?>.  The  bulk  relaxation  mechanisms 
only  affect  the  magnitude  of  <t»§>,  not  its  spatial  distribution. 

Minguzzi  et  al.3®  have  shown  that  the  spatial  distribution  of  <!•$>  is, 
in  general,  characterized  by  a  sum  of  diffusion  modes.  In  the  case  of  the  gas 


Hi 
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cell  clock,  though,  where  the  optical  pumping  rate  is  low  and  the  resonance 
cell  is  fairly  uniformly  illuminated,  a  reasonable  approximation  is  to 
consider  only  the  first  order  mode.  Thus,  we  have  for  the  axial  distribution 
of  <t*£> 

<f»£>  (z)  *  <!•?>  sin  (irz/i)  *  (19) 

P 


Because  is  simply  related  to  <!•$>  by  the  equation^ 


1  +  1 

2  ‘  g 


(20) 


we  have  for  the  spatial  distribution  of  n  *  considering  only  the  first  order 

A 

axial  diffusion  mode. 


Vz) 


H  -  r  > *in  <r> 


(21) 


In  essence,  we  have  decomposed  the  contributions  of  diffusion  to  n.  into  two 

A 

parts.  In  the  first  part,  we  assume  that  diffusion  lias  some  limiting  effect 

on  the  magnitude  of  n  .  and  that  this  effect  can  be  Incorporated  into  the  gen- 

A 

eralized  Vanier  theory  through  contributions  to  y  and  y  .  In  the  second 

«  1  z 

part,  we  recognize  that  the  final  form  of  the  axial  distribution  of  n  is  to 

A 

some  degree  known,  so  a  functional  form  can  be  superimposed  onto  the  solution 
n^.  In  this  way,  the  major  physical  aspects  of  diffusion  are  brought  into  the 
model  without  recourse  to  cumbersome  and  obscure  expressions. 

The  bulk  relaxation  mechanisms  which  contribute  to  y^  and  y ^  are 
typically  buffer  gas  and  spin  exchange  collisions.  Additionally,  for  com¬ 
pleteness,  one  should  Incorporate  the  slowdown  factor,  (61  +  1 ) / ( 81  +  4),  into 

25 

the  spin  exchange  contribution  to  the  transverse  relaxation  rate  y^.  Thus, 

if  we  define  a__  and  a  as  the  buffer  gas  relaxation  and  spin  exchange  cross 
do  ex 
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a 


sections,  respectively, the  phenomenological  relaxation  rates  y^  and  y. 


become 


D/k  +  [BG]  va_„  +  [Rb]va 

OO  GX 


(22a) 


D/R  +  [BG]vo, 


+  (4rH)  [Rblvot 


(22b) 


where  [BG]  is  the  number  density  of  buffer  gas  molecules  or  atoms,  and  v  is 

the  average  relative  velocity  between  colliding  species. 

Finally,  by  assuming  that  the  Rb  atoms  are  essentially  frozen  in  place 

because  of  the  presence  of  the  buffer  gas,  we  can  determine  n'(z)  by 

A 

substitution  of  the  position  dependent  photon  absorption  rates,  A  and  B,  along 
with  the  suitable  expression  describing  u>j(z)  into  Eq.  (15).  For  typical  gas 
cell  frequency  standards  employing  a  TEqj j  or  TEjjj  microwave  cavity  this 
expression  is  simply 


Jl(z)  "  “lp  8in  (?) 


where  w  is  the  peak  microwave  Rabi  frequency  in  the  cavity. 

4.  DETERMINATION  OF  THE  DEPENDENCE  OF  PHOTON  ABSORPTION  RATES  ON 
AXIAL  POSITION 

To  calculate  the  positional  dependencies  of  the  optical  photon  absorption 
rates,  we  first  assume  that  the  axial  optical  intensity  distribution  can  be 
approximated  by  a  Bouger-Lambert  law  of  the  form 


SV-Y-  Y-'lv" 


(24a) 


IjU.z)  =  IjQ  exp  (-{(1  -  c)  [Rb87]  +  (|j)  [Rb85]}ajZ)  (j  *  3,  4)  (24b) 


where  C  Is  to  be  thought  of  as  a  resonance  cell  parameter  describing  the 

■'global"  fractional  population  in  the  F  *  1  hyperfine  multiplet  as  a  result 

of  optical  pumping.  In  effect,  5  is  a  zeroth  order  approximation  to 

(z).  Feu:  the  low  optical  pumping  rates  expected  from  lamps,  this 

approximation  should  be  a  reasonable  description  of  reality,  because  in  the 

limit  of  negligible  optical  pumping,  rT  (z)  in  Eqs.  (12a)  and  (12b)  is 

A 

independent  of  the  light  intensity,  and  the  standard  Bouger-Lambert  law  is 
recovered.  We  note,  however,  that  recent  theoretical  work  has  indicated  that 
this  approximation  may  not  be  generally  valid. 22,23  Equations  (24a)  and  (24b) 
will  permit  us  to  calculate  the  positionally  dependent  optical  photon 
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1  I.(C,  *)o.  (25a) 

j-3  J  3 

2 

l  I . (?  ,  z)a  (25b) 
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absorption  rates  as  functions  of  f. 


A(C.  *>  *  <fc> 
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The  task  which  now  needs  to  be  discussed  is  the  determination  of  c . 
Considering  Eqs.  (14)  and  (24),  it  is  clear  that  any  consistent  choice  of  5 
must  satisfy  the  equation 


I  MS,*)  *  l  exP  {-P„  tRb87]o  l) 

k-1  i-1  10 
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(26) 


The  simplest  choice  of  5  which  satisfies  this  requirement  is 


n  *  -i 

®  1 


7  / 
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Substituting  Eqs.  (15)  and  (21)  into  Eq.  (27),  our  expression  for  5  becomes 


gK  SkS,  1  "  A^C)I  sin  (f5-) 

c  -F  +  — /  {  gjBU)  +  Yij  “gbiA(5r+77i}  dz "  ° 


where  we  have  explicitly  noted  that  A  and  B  depend  parametrically  on  5.  Using 
a  Newton-Raphson  method, ^  we  can  solve  Eq.  (28)  numerically  for  £.  A 
discussion  of  the  results  of  this  procedure  for  determining  5  and  their 
consistency  will  be  presented  subsequently. 

In  summary,  with  the  determination  of  5,  A(z)  and  B(z)  are  specified  and 

when  substituted  into  Eq.  (15),  determine  n!»  The  terms  n  and  TT  follow 

A  A  A 

directly  from  n'  via  Eqs.  (21)  and  (13),  respectively.  Substituting  n  into 
A  A 

Eq.  (14)  yields  the  transmitted  light  intensity  I(A ,  l  )  and  Che  optically 
detected  hyperfine  lineshape,  S(A)  ■  I(»,  t  )  -  I(A»  l  ). 

C.  INTEGRATION  OF  ATOMIC  PHYSICS  WITH  FEEDBACK  ANALYSIS 

It  is  a  straightforward  process  to  extract  the  quantities  required  by  the 
feedback  analysis  from  the  preceding  study  of  the  physics  occurring  in  the 
absorption  cell.  If  <  and  A^  are  the  clock  photodiode's  sensitivity  (A/W)  and 
area,  which  is  assumed  to  be  identical  Co  Che  clock  cavity's  cross-sectional 
area,  the  total  lamp  induced  photocurrent  is  given  by 


i  ■  <A  I(",  O 
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The  slope  of  the  discrimination  pattern  is  given  by 

m  -  5  1“(0,  1)  (30) 

in 

where  26^  is  the  modulation  depth  of  the  microwave  frequency,  and  I"(0,  i)  is 

the  second  derivative  of  the  transmitted  light  intensity  evaluated  at  the 

center  frequency  of  the  hyperfine  resonance.  By  equating  the  slope  of  the 

discriminator  pattern  with  5  I”(0,  l) ,  we  have  made  the  small-modulation- 

m 

depth/quasi-static  approximation  with  regards  to  the  effect  of  the  microwave 
frequency  modulation  on  the  signal  lineshape.^  As  long  as  we  confine  our¬ 
selves  to  frequency  modulation  regimes  such  that  5  /!*_  ~1  and  u  /r,  ~1 , 

where  a»  is  the  modulation  frequency,  then  the  recent  work  of  Audoln  et  al.  9 
M 

would  indicate  that  this  approximation  is  reasonably  valid.  As  will  be  demon¬ 
strated  in  the  following  paragraphs,  the  exclusion  of  higher  order  modulation 
effects  does  not  appear  to  have  affected  the  ability  of  our  model  to  accurate¬ 
ly  predict  the  performance  of  current  generation  Rb  clocks.  It  is  not  clear, 
though,  that  this  would  be  true  under  all  potential  clock  operating  condi¬ 
tions.  If  clock  conditions  varied  significantly  from  nominal  values  used  in 
current  standards,  the  ramifications  of  microwave  modulation  would  have  to  be 
reconsidered.  However,  these  effects  could  be  approximated  in  the  present 
model  by  using  the  work  of  Audoln  et  al.  to  determine  a  correction  factor 

for  Eq.  (30).  In  our  calculations,  6  is  considered  equal  to  the  frequency 

in 

separation  of  the  inflection  points  of  the  approximately  Lorentzlan  resonance 

lineshape,  and  it  is  assumed  that  u  <  5  . 
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III.  MODEL  VERIFICATION 


In  the  model  of  the  clock  signal  lineshape  just  presented,  quite  a  few 

assumptions  were  made  in  order  to  make  the  problem  tractable.  Several  of 

these  assumptions  are  Inherent  to  the  generalized  Vanier  theory  and  have  been 

on 

validated  by  recent  experimentation.  These  can  be  considered  as  microscopic 

assumptions,  because  they  are  related  to  the  statistical  behavior  of 

individual  atoms.  The  other  assumptions  in  the  model  (i.e.,  the  functional 

form  of  the  axial  light  intensity  distribution,  and  the  axial  distribution 

of  n  as  a  result  of  diffusion)  should  be  considered  as  macroscopic 
A 

assumptions,  because  they  are  related  to  the  alkali  vapor  as  a  whole.  We  are 
on  less  firm  ground  experimentally  with  these  assumptions.  Thus,  before 
proceeding  to  calculations  of  clock  performance,  it  would  be  prudent  to 
explore  the  validity  of  these  assumptions.  This  will  be  accomplished  by 
showing  that  the  model,  in  its  present  form,  leads  to  qualitatively  correct 
predictions  of  the  global  fractional  population  g,  the  axial  distribution 
of  <!•?>,  and  the  signal  lineshape  as  functions  of  temperature. 

The  global  fractional  population  g  is  actually  a  quite  useful  parameter 
for  understanding  gas  cell  clock  performance,  because  it  reflects  the  degree 
of  optical  pumping,  and  hence  the  clock  signal  amplitude,  for  the  resonance 
cell  as  a  whole.  Qualitatively,  one  expects  g  to  depend  on  temperature  in  two 
ways:  (1)  as  the  resonance  cell  temperature  rises,  the  global  optical  pumping 
rate  necessarily  changes;  and  (2)  the  longitudinal  relaxation  rate  is  tempera¬ 
ture  dependent  as  a  result  of  spin  exchange.  Thus,  for  high  temperatures,  one 
expects  g  to  approach  its  statistical  value  g^/g.  At  low  temperatures,  when 
the  vapor  is  optically  thin,  one  should  be  able  to  arrive  at  a  relatively 
simple  expression  for  g  because  in  this  case,  one  can  assume  that  A  and  B  are 
constants  independent  of  axial  position.  Regarding  Eq.  (28)  in  this  optically 
thin  regime,  one  has 


t  8b  u  VB  ~ A) _ 

thin  g  '  w[gfl(B  +  Yj )  +  gb(A  +  Y^r 
which  is  still  a  temperature  dependent  quantity  because  of  y 
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If  we  now  consider  Che  case  of  optical  puraping  out  of  only  Che  F  -  b 

hyperflne  raultiplet,  as  a  specific  example,  Che  preceding  arguments  would 

require  us  Co  make  Che  following  predictions  regarding  (I)  at  low  resonance 

cell  ceraperacures,  c  should  be  relatively  small;  (2)  at  high  temperatures,  z 

should  saturate  at  3/8  in  the  case  of  Rb;  and  (3)  for  intermediate 

temperatures  one  should  have  z  >  t  ,  .  ,  because  r  .  increases  as  a 

correct  thin*  *thin 

function  of  temperature  as  a  result  of  only  one  mechanism,  whereas  t 

correct 

increases  as  a  result  of  two  cooperating  mechanisms.  A  sample  calculation 
of  c  over  the  temperature  range  10-120*C  is  presented  in  Fig.  3.  The 
qualitative  agreement  between  our  expectations  and  the  calculations  lends 
confidence  to  the  Interpretation  and  use  of  z  in  the  clock  signal  model. 

Because  the  atoms  are  essentially  frozen  in  place  as  a  result  of  the 
presence  of  the  buffer  gas,  the  maximum  change  In  light  intensity  per  unit 
length  (for  microwaves  on  and  off  resonance)  occurs  in  a  unique  spatial  region 
of  the  resonance  cell  determined  by  the  light  intensity  distribution,  the 
microwave  field  distribution,  and  diffusion.  The  ability  to  determine  this 
unique  spatial  region  is  critically  important  for  a  model  of  gas  cell  clock 
performance,  because  the  microwave  field  strength  and  light  intensity  in  this 
region  have  the  primary  Influence  on  both  the  clock's  Q  and  signal-to-noise 
ratio.  Furthermore,  variations  in  this  position  with  microwave  power  and 
light  intensity  are  responsible  for  the  so-called  position  shift^z  and 
inhomogeneous  light  shift, respectively. 

As  previously  discussed,  the  present  model  is  one  dimensional.  Thus,  if 
the  model  is  to  be  trusted  in  clock  performance  predictions,  it  should  be  able 
to  predict  the  qualitatively  correct  axial  distribution  of  <1»$  >  as  a 
function  of  resonance  cell  temperature.  For  the  case  of  an  optically  thin 
cell  and  relatively  low  optical  pumping  light  intensities,  one  would  expect  a 
distribution  of  >  peaked  in  the  center  of  the  resonance  cell.  As  the 

resonance  cell  temperature  was  raised,  and  as  the  cell  became  more  opaque  to 
the  optical  pumping  light,  one  would  expect  the  position  corresponding  to  the 
maximum  <!•§  >  to  shift  towards  the  front  of  the  resonance  cell.  In  Fig.  4, 
we  have  plotted  the  model's  predictions  of  <!•§  >  as  a  function  of  axial 
position  using  the  parameters  listed  in  Table  1.  We  note  that  the  curves  are 


GLOBAL"  FRACTIONAL  POPULATION  IN 
LOWER  HYPERFINE  MULTIPLES  f 


RESONANCE  CELL  TEMPERATURE  (°CI 


Fig.  3.  Global  Fractional  Population  ?  in  the  F  ■  b  *  I  -  1/2  Hyperflne 

Multiplet  as  a  Function  of  Resonance  Cell  Temperature  for  the  Case  of 
Optical  Pumping  Out  of  Only  This  Multiplet.  The  solid  line  is  a  thin 
vapor  approximation  where  the  spin  exchange  rate  is  the  only 
temperature  dependent  quantity;  dots  correspond  to  the  actual  case 
where  the  optical  thickness  of  the  vapor  must  also  be  considered.  In 
the  limit  of  no  optical  pumping,  the  equilibrium  value  of  c  Is  3/8. 
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The  Measure  of  Hyperfine  Polarization  <!•§>  has  an  Axial  Spatial 
Distribution  That  is  Temperature  Dependent.  For  relatively  low 
resonance  cell  temperatures,  the  spatial  distribution  is  primarily 
determined  by  diffusion.  As  the  resonance  cell's  alkali-metal  vapor 
becomes  progressively  more  optically  thick,  efficient  optical  pumping 
becomes  confined  to  regions  closer  to  the  entrance  window  of  the 
resonance  cell. 


•«*j 


Table  1.  Parameters  for  Clock  Signal  Model 


Parameter 

Value(a) 

Uncertainty 

Linewidth  of 

2  GHz 

±0.2  GHz 

lamp  lines 

Temperature 

65°C 

±  5°C 

Photodiode 

responsivity 

0.5  A/W 

Ratio  of  buffer 

0.64 

gas  light  to 

Rb  resonance  light 

P0(DpF  -  1) 

0.029  mW 

±  5X 

po(d2,f  -  1) 

0.046  mW 

P0(DpF  -  2) 

0.018  mW 

*5%  (b) 

Pq(D2,F  -  2) 

0.029  mW 

Peak  Rabi 

450  Hz 

frequency  u>^ 

Buffer  gas 

n2 

Buffer  gas 

10  Torr 

pressure 

Light  shift 

-3  x  10"10  (c) 

±  1  x  10"10 

coefficient 

Lamp  Intensity  1.2  x 

I0'9  +  (c) 

spectral  density 

f 

(a)Ref.  9 


(b)  Optical  powers  were  varied  in  the  calculations  according  to  the  rule  that 
if  p0(Di»F  *  1)  changed  by  +  5Z,  then  simultaneously  PQ(Di*P  *  2)  changed 
by  *  5Z. 


(c) 


Ref.  16 
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very  similar  to  those  of  Minguzzi  et  al.,  where  the  effects  of  diffusion 
were  analyzed  more  rigorously,  and  are  again  in  qualitative  agreement  with  our 
expectations. 

Of  course,  the  most  critical  test  for  our  macroscopic  assumptions 
concerns  the  signal  lineshape.  Experimentally,  it  is  known  that  the 
transmission  signal  resonance  is  very  well  approximated  by  a  Lorentzian, 
being,  however,  somewhat  more  sharply  peaked. ^  Regarding  the  model,  there  is 
no  reason  a  priori  to  expect  a  lineshape  resembling  anything  like  a 
Lorentzian;  yet  if  this  is  not  predicted  over  a  reasonable  range  of  operating 
temperatures,  the  validity  of  the  model  must  be  questioned.  In  Figs.  5a  and 
5b,  we  show  the  signal  lineshape  calculated  at  20  and  100#C,  respectively, 
using  the  parameters  listed  in  Table  1.  The  dots  are  model  calculations,  and 
the  solid  curve  is  a  least  squares  fit  of  the  model  results  to  a  Lorentzian. 

As  is  readily  apparent,  the  model  gives  excellent  agreement  in  the  wings  of 
the  lineshape,  but  near  resonance  is  more  sharply  peaked,  as  expected.  These 
predictions  are  in  qualitative  agreement  with  the  experimental  results  of 
Camparo  and  Frueholz.^  Thus,  our  macroscopic  assumptions  appear  to  be  valid, 
at  least  over  normal  clock  operating  conditions,  permitting  the  final  test  of 
the  model:  a  comparison  between  predicted  and  measured  Rb  gas  cell  clock 
performance. 

9  10 

Recently,  Riley  and  Lynch  and  Riley  performed  a  stability  test  on  a 
newly  designed  Rb  frequency  standard.  This  test  was  somewhat  unique,  because 
many  of  the  input  parameters  required  by  the  present  model  of  clock  perform¬ 
ance  were  documented.  Thus,  these  results  can  be  used  to  establish  the 
accuracy  of  the  present  model  in  predicting  clock  performance.  The  parameters 
required  for  the  calculation  are  listed  in  Table  1,  along  with  reasonable 
estimates  of  their  values,  based  on  the  information  presented  by  Riley  or  on 
our  own  measurements  on  similar  components  of  the  Rb  clock.  The  uncertainties 
presented  with  some  of  the  more  critical  parameters  represent  our  best  esti¬ 
mates  of  the  inaccuracies  associated  with  extracting  the  parameters  from  the 
information  presented  by  Riley.  In  particular,  the  uncertainty  associated 
with  the  resonance  cell  temperature  represents  the  difficulty  associated  with 
measuring  liquid  alkali  metal  temperatures  accurately,  and  the  difficulty  of 
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OpticBUy  Detected  Hyperfine  Transition  Lineshape.  (a)  T  -  20°C, 

(b)  T  -  lOO-C.  The  dots  are  the  calculated  change  in  transmitted 
light  intensity  as  a  function  of  microwave  frequency;  the  solid  curve 
is  a  least  squares  Lorentzian  fit.  Note  that  at  reasonable  resonance 
cell  temperatures,  the  lineshapes  are  well  approximated  by  a 
Lorentzian;  they  are,  however,  somewhat  more  sharply  peaked. 


accurately  converting  these  liquid  temperatures  to  saturated  alkali  vapor 
pressures.  The  only  parameter  that  could  not  be  determined  by  the  information 
presented  by  Riley  was  the  peak  Rabi  frequency.  This  parameter  was, 
therefore,  evaluated  by  matching  the  calculated  signal  linewidth  with  Riley's 
measured  signal  linewidth.  We  note,  however,  that  this  parameter  can  be 
determined  experimentally  quite  easily  by  the  adiabatic  rapid  passage  method 
of  Frueholz  and  Caraparo.  The  choice  of  buffer  gas  was  not  crucial,  so,  for 
convenience,  nitrogen  at  10  Torr  was  used. 

Table  2  and  Fig.  6  display  the  results  of  our  calculations  along  with 
uncertainties  in  the  calculated  quantities  which  are  associated  with  the 
uncertainties  in  the  input  parameters.  Several  points  concerning  the  pre¬ 
dicted  results  are  worth  noting.  First,  the  predicted  Q  and  discriminator 
slope  agree  very  well  with  the  measured  values,  especially  because  the 
measured  discriminator  slope  is  probably  a  bit  too  large.  Based  on  the 
measured  Allan  variance  and  DC  photocurrent,  one  would  have  expected  to 
measure  a  discriminator  slope  of  268  pA/Hz.  Furthermore,  the  results  indicate 
that  the  line  Q  is  relatively  insensitive  to  reasonable  variations  in  the 
global  optical  pumping  rate,  but  that  the  signal  amplitude  is  very  sensitive 
to  these  changes.  This  is  precisely  the  same  conclusion  reached  by  Matsuda  et 
al.  in  their  experimental  studies  of  the  signal  characteristics  of  the  Rb  gas 
cell  clock. ^  Finally,  the  calculated  Allan  variance  agrees  very  well  with 
that  measured  by  Lynch  and  Riley, as  shown  by  Fig.  5.  Thus,  within  the 
uncertainties  associated  with  the  input  parameters,  we  have  demonstrated  the 
validity  of  a  model  of  gas  cell  clock  performance  that  is  completely  independ¬ 
ent  of  measurements  of  the  clock  signal's  Q  or  signal-to-noise  ratio. 

Vanler  et  al.  have  carried  out  extensive  analyses  of  the  effect  of  loop 
attack  time  on  clock  performance.  As  an  additional  demonstration  of  the  model 
capabilities,  the  Allan  variances  for  two  different  servo-control  loop  attack 
times  are  calculated.  The  crystal  oscillator  is  assumed  to  have  a  power 
spectral  density  of  fractional  frequency  fluctuations  of*2 


(f) 


6  x  10 


-26 


+  1  x 


10"26f  +  2  x 


io-28f2 
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Table  2.  Comparison  of  Predicted  Clock  Performance 
with  Measured  Clock  Performance 


Quantity 


Predicted  Valve  Measured  Valve 


Fractional  Population  ; 

0.359 

Fraction  Population 

Change:  n  -5 

0 

4.3  x  10" 

•3 

DC  photocurrent 

81.8  uA 

82  uA 

Photocurrent  change: 
idc  '  1a*0 

+352 

(45.7  ) 

-302 

nA 

Q: 

+1.22 

(1.63  ) 

-1.82 

x  107 

1.6  x  107 

Discriminator  pattern 

slope:  6  I"(0,O 
in 

+93 

(251  ) 

-74 

pA/Hz 

351  pA/Hz 

Short-Term  Allan 
variance:  0  (t) 

V 

+1.1 

(3.1  )x 

-0.9 

10"12/  /7 

2.8  x  10-12 

LOG  [AVERAGING  TIME  t  (sec)] 


Fig.  6 


.  Frequency  Stability  of  the  EG&G  Rb  Gas  Cell  Frequency  Standard  In 
Terms  of  Allan  Variance  a  (t).  The  dashed  line  is  the  standard's 
experimental  performance ,^while  the  solid  line  is  the  predicted 
performance  yielded  by  the  model.  The  error  bar  represents  the 
uncertainty  in  the  standard's  predicted  performance  due  to 
uncertainties  in  the  input  parameters.  For  averaging  times  less  than 
104sec,  the  standard's  performance  is  limited  by  shot  noise  at  the 
photodetector.  When  the  averaging  time  is  greater  than  104sec,  clock 
performance  is  limited  by  discharge  lamp  intensity  fluctuations 
transferred  into  frequency  fluctuations  through  the  light  shift 
effect. 


while  the  atomic  system  behaves  as  previously  calculated.  Two  attack  times, 
0.1  and  10  sec,  are  considered.  The  resultant  Allan  variances  are  displayed 
in  Fig.  7.  These  results  were  calculated  quite  efficiently  because  of  the 
analytical  evaluation  of  the  integrals  of  Eq.  (8).  As  an  aside,  we  note  that 
the  present  analysis  of  the  hyperfine  lineshape,  based  on  the  quasi-static  ap¬ 
proximation,  will  not  be  valid  at  extremely  short  attack  times.  At  very  short 
attack  times  the  Rb  atoms  will  no  longer  be  able  to  respond  adiabatically  to 
the  microwave  frequency  variation  induced  by  the  control  loop.  Furthermore, 
at  short  attack  times,  atomic  coherence  effects  should  be  considered.  For 
typical  gas  cell  frequency  standards,  a  reasonable  limit  to  the  validity  of 
our  model  would  be  an  attack  time  of  0.01  sec. 


-2  -1  0  1  2  3  4  5 

LOG  [AVERAGING  TIME  r  (sec)] 

87 

ated  Allan  Variance  a  (t)  for  Passive  Rb  Gas  Cell  Standard 

y 

ro  Control  Loop  Attack'Times .  Crystal  oscillator  performance  as 
:ied  in  Ref.  17.  The  atomic  system  performance  is  that  of  the 
>rototype  standard. 


IV.  CONCLUSIONS  AND  SUMMARY 


In  practice,  a  clock  model's  usefulness  is  intimately  connected  to  the 
realization  of  a  device  with  optimal  performance  characteristics  for  a 
specific  application.  Thus,  though  reasonable  accuracy  in  model  predictions 
is  required,  equally  Important  attributes  of  a  clock  model  are  its  ability 
to:  (I)  analyze  design  trade-offs  for  a  specific  application,  (2)  diagnose 
the  clock  parameters  that  limit  performance,  and  (3)  predict  the  avenues  for 
further  improvement.  Because  the  present  model  is  inherently  numerically 
faster  than  its  predecessors,  as  a  result  of  the  analytic  evaluation  of  the 
spectral  density  integrals  and  density  matrix  Inversion,  design  trade-offs  can 
be  quickly  evaluated.  Furthermore,  because  the  model  treats  the  physics  of 
the  resonance  cell  signal  in  a  semirlgorous  yet  perspicuous  manner,  diagnoses 
of  performance  limits  and  predictions  for  improvements  are  in  terms  of 
physically  meaningful  and  conceptually  lucid  parameters. 

As  an  example  of  the  present  model's  ability  to  predict  possible  avenues 
of  improvement  for  the  Rb  gas  cell  clock,  we  have  computed  the  short-term 
performance  for  a  clock  with  the  parameters  listed  in  Table  1  as  a  function  of 
microwave  power,  or  what  is  equivalent  to  the  peak  Rabl  frequency  in  the  clock 
cavity.  These  results  are  shown  in  Fig.  8.  As  is  clear  from  the  figure,  the 
model  predicts  that  a  short-term  stability  of  6  x  10-l^//x  could  be  achieved 
by  reducing  the  6835  MHz  microwave  power  fed  into  the  cavity  by  23  dB.  Of 
course,  these  predictions  must  be  tempered  by  the  realization  that,  for  these 
Rabi  frequencies,  one  might  not  be  in  the  small-modulation-depth/quasi-state 
regime  with  regard  to  the  discriminator  slope,  and  that  the  microwave  field 
strength  displays  only  a  longitudinal  variation  within  the  absorption  cell. 
However,  the  predictions  do  indicate  a  possible  avenue  for  improvement  that 
must  be  experimentally  verified. 

Finally,  with  a  straightforward  modification,  the  microwave  detuning  in 
Eq.  (9)  can  be  made  to  depend  upon  local  magnetic  field  strengths,  local 
optical  pumping  light  intensities,  the  buffer  gas  pressure,  the  alkali  vapor 
pressure,  etc.  The  model  lineshape  would  then  exhibit  the  Inhomogeneous 
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Fig.  8.  Short-Term  Allan  Variance  Slope  for  the  Parameters  of  Table  1  as  a 
Function  of  Microwave  Rabl  Frequency,  or  Equivalently,  Microwave 
Power  Fed  into  the  Cavity  Assuming  a  Q  ■  100.  These  calculations 
would  indicate  roughly  a  factor  of  5  improvement  in  the  short-term 
performance  by  decreasing  the  microwave  power  by  ~  23  dB. 
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nature  that  results  in  the  position  shift  effect^  an<j  t^e  inhomogeneous  light 
shift. ^  With  the  Inclusion  of  these  resonant  frequency  shift  terras,  the 
model  could  be  used  to  Investigate  possible  mechanisms  of  gas  cell  dock, 
frequency  drift,  which  is  a  principal  impediment  to  the  gas  cell  frequency 
standard's  more  widespread  use. 
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APPENDIX:  RESULTS  OF  COMPLEX  CONTOUR  INTEGRATION 


When  extracting  the  Allan  variance  frora  the  spectral  density  of  frac¬ 
tional  frequency  fluctuations,  integrals  of  the  form 


A(t)  ■  /  - - - t  sin*  (irrf)  df 

0  1  +  (f/f  ) 

n 

with 


(Al) 


a  -  -4,  -2,  -I,  1,  2 


must  be  evaluated.  Evaluation  of  these  Integrals  may  be  performed  using 
complex  contour  techniques.  To  facilitate  evaluation,  sin^(yf)  was  written  in 
terms  of  complex  exponentials.  For  the  even  integrands,  the  contour  selected 
extended  from  -R  to  +R  along  the  real  axis,  with  a  semicircular  Indention 
around  origin  when  required.  The  contour  was  closed  along  the  upper  half  of 
the  circle  |z|  »  R  in  the  complex  plane  with  R  ultimately  permitted  to  go  to 
infinity.  For  odd  integrands,  the  contour  was  composed  of  three  segments. 

The  first  segment  extended  from  MR  to  zero  along  the  Imaginary  axis,  with 
semicircular  Indentions  about  singularities.  The  second  segment  extended  from 
0  to  +R  along  the  real  axis,  followed  by  the  third  segment,  a  quarter  of  the 
circle  |z|  ■  R  connecting  the  other  segments.  Again  R  was  permitted  to  go  to 
infinity.  In  evaluating  the  integral  for  o  •  -  1,  it  was  found  to  be 
expeditious  to  first  differentiate  tne  integrand  with  respect  to  y(“xt), 
eliminating  the  singularity  at  the  origin.  The  resultant  integral  was 
evaluated.  Indefinite  integration  with  respect  to  y  was  then  performed. 

The  results  of  this  analysis  are  as  follows: 

,  ,  .  ,  '4"fn  -2„f 

a  -  -  4:  A(t)  -  j  [y  (»T>*  -  J  (| - €  n  +  3/4)]  (A2) 

^n  ^n 


51 


,N 


X^s'CvT'wT-.Y-o 


J*. 


tAiXw^ 


2  -4irrf  -2irtf 

a  -  -  2:  A(t)  -  ~  (e  n  -  4e  n  +  3) 

n 


a  -  -  1:  A(t)  -  -fir  [F(4f n*T)  -  2F(2fnirr)] 

n  1 


(A4a) 


F(cirt) 


[EK-cirr)  - 


n  /  xn-k 

n  r  (-1)  (cirt)  i 


<ct)  *!  (-.)"  J  ; 


Bn,. n-k 
c  r  (  city  ) 


+  l -  [Ei(cwr)  -  e  EI(-cxt) «+  \  \  k),1 

n-1  k-1  ^ 

,  A  B 

-  -  (  £  i  +  l  -) 

c  L.  n  L,  n' 


(A4b) 


for  cwt  <  20.  Ei(x)  is  the  exponential  Integral  function,  and  A  and  B  are 
selected  to  ensure  convergence  of  the  respective  sums.  When  cirt  >  20,  various 
asymptotic  expansions  may  be  employed  In  deriving  F(cirx)  to  yield 


M  (k  -  2)! 

F(cxt)  -  -=•  [ln( cirr)-  l  - r— pi 

c  k-3  (cirt)* 

k  odd 


(A4c) 


with  M  selected  to  ensure  series  convergence  to  desired  level  of  accuracy. 

When  o  equals  either  +1  or  +2,  the  integral  of  Eq.  (Al)  is  not  finite. 
One  way  this  can  be  overcome  is  by  limiting  the  range  of  integration  to  a 
raaxmium  frequency.  •  Physically  this  represents  the  maximum  frequency  to  which 
the  measurement  system  responds.  An  alternative  means  of  addressing  this 


Handbook  of  Mathematical  Functions,  M«  Abramowitz  and  I.  A.  Stegun,  eds., 
Dover  Publications,  Inc.  (1965),  pp.  227-251. 
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problem,  preserving  analytic  evaluation  of  the  integrals,  is  to  modify  the 
integrand  in  such  a  manner  that  it  goes  smoothly  to  zero  above  a  cutoff 
frequency  fc>  This  is  carried  out  by  altering  Eq.  (Al)  to  give 


A(x)  *  /  - ^ - j  - - - j  sin^(irTf)df 

0  1  +  <f/f  )  1  +  (f/f  ) 

n  c 


Evaluating  Eq.  AS  for  a  ■  +1  and  +2  yields 


a  -  +1:  A(t)  -  - — - = —  {G(4irf  t)  -  G(4irf  r)  -  4[G(2irf  t)  -  G(2irf  t)] 

16(f  2rf  2)  C  n  C  .  n 

c  n 


+  6  in(f  / f  )} 
c  n 


with  G(x)  -  e“x  Ei(x)  +  ex  Ei(-x) 


a  ■  +2:  A(t) 


w(f  )  -4wfT 

— — — r1£c(e 


-2irf  t 

-4e  c  +3) 


-4irf  t  -2irf_T 


-fn  (e  “  -4e 


+  3)] 


LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for 
national  security  projects,  specializing  In  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  Investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  Investigations  la  the  technical  staff's  wide-ranging  expertise  and  Its 
ability  to  stay  current  with  new  developments.  This  expertise  Is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  Tluld  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  exclmer  laser  development  Including  chemical  klnetlca, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  statm-apeclf lc  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of -f ield-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemiatry ,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Prograu  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  Intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  mlcrowave/mllllmeter  wave  measurements, 
diagnostics  and  radlometry,  mlcrowave/mllllmeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  In  space  and  enemy-induced 
environments. 

Space  Sclerices  Laboratory:  Magnetospherlc,  auroral  and  cosmic  ray 
physics,  wave-particle  Interactions,  magnetospherlc  plasma  waves;  atmospheric 
and  Ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  Infrared  astronomy. 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  Ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
Instrumentation. 


